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Abstract 

f — We consider the generalised root identities introduced in [1] for simple 

functions, and also for r(z + 1) and (^(s). Where these were almost ex- 
^ ^ ' clusively examined just for /i G Z in [1], here we focus exclusively on the 

L case of jj,^ Z. For the simplest function f{z) = z, and hence for arbitrary 

polynomials, we show that they are satisfied for arbitrary /x 6 R (and 
^ hence for arbitrary /i G C by analytic continuation). Using this, we then 

develop an asymptotic formula for the derivative side of the root identities 
2 for r{z + 1) at arbitrary /x G R, from which we are able to demonstrate 

C numerically that r{z + 1) also satisfies the generalised root identities for 

' ' arbitrary n, not just integer values as verified in [1]. Finally, using the 

working in [1] we examine the generalised root idcntitcs for also for non- 
integer values of /i. Having shown in [1] that ( satisfies these identities 
exactly for jj, G Z>2 (and also for // = 1 after removal of an obstruction), 
in this paper we present strong numerical evidence first that ( satisfies 
them for arbitrary fi > 1 where the root side is classically convergent, 
and then that this continues to be true also for — 1 < /i < 1 where Ce- 
saro divergences must be removed and Cesaro averaging of the residual 
partial-sum functions is required (when fx < 0). A careful examination 
of the calculations in a neighbourhood of = also sheds light on the 
appearance of the 2d In-divergence that was handled heuristically in [1] 
and why the assignment of 2d Cesaro limit to this in [1] is justified. The 
numerical calculations for /i > are bundled in portable R-code which 
can be readily used to further verify these claims; the code for the case 
f I —1 < /i < 0, including the Cesaro averaging required when /i < 0, is in 

VBA in an XL spreadsheet which can also be used to check these claims 
and conduct further tests. Both the R-scripts and XL spreadsheet are 
made available with this paper, along with supporting files. 
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1 Introduction 

In [1] we introduced the generalised root identities, namely: 



r(/u) \dz 



{zq— roots ri} 
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The LHS in (1) is the derivative side of the root identities for / at zq and /i; the 
RHS is the root side. Writing (i(so,/i) for the expression on the derivative side 
and r{so,^l) for the root side, saying that a function / satisfies the generalised 
root identities means that (i(so, /i) and r(,so, /i) are identically equal as functions 
of two complex variables (we may omit reference to the underlying function / 
in the notation for d and r unless the context requires it). 

In general, within an equivalence class of functions (equivalent up to mul- 
tiplication by a nowhere-zero entire function), we expect there to be a unique 
representative which satisfies these root identities. For many functions of great 
interest, however, we saw in [1] that they themselves satisfy the generalised root 
identities, at least for many values of sq and /U, without any need to go exploring 
within their equivalence classes. For example in [1] we showed that: 

(a) the function f{z) = cos(^) satisfies the generalised root identities at 
Zq = for arbitrary /x e C, with the root identities in this case being precisely 
equivalent to the famous functional equation for C in /x. 

(b) the function r(z + l) satisfies the generalised root identities for arbitrary 
Zo € C and /x G Z, with a renormalisation being required when /x = 1, and with 
a distributional interpretation of the derivative side and a Cesaro evaluation of 
the root side required when /i G Z<o. 

(c) the function ({s) satisfies the generalised root identities for arbitrary 
So S C and e Z>2, while n~i({s) satisfies the identities when = 1 after the 
removal of an "obstruction" in this instance (the resulting /i = 1 identity being 
equivalent to the Hadamard product formula for (). Additionally we showed 
that the derivative side of (1) should be identically zero as a function of sq when 
H G 1<Xi and we conducted three Cesaro calculations of the root side for such 
cases, namely for /x = 0, /i = — 1 and \i = —2. For = and /U = — 1 we 
verified the claimed result, at least modulo an estimate on the argument of the 
zeta function, S(T), in the case /i = — 1 which is known to hold if the Riemann 
hypothesis (RH) is true. For fj, = —2, however, we claim that the resulting 
Cesaro evaluation of the root side suggests that the Riemann hypothesis must 
be false. 

As discussed in [1], however, the fact that these computations were conducted 
for isolated values of n, rather than as a systematic extension in open regions 
of the /u-plane, is a deficiency of the argumentation thus far (in [1]. the case 
of f{z) = cos(^) when = is in fact the only instance in the whole paper 
where the case of arbitrary /i ^ Z is considered for root identity calculations) . 

In this paper, therefore, we address this by explicitly focusing on the root 
identities for arbitrary n and zq in M (and then C by analytic continuation), 
verifying them in this general case also for the functions considered in [1], in 
particular F and (. 

As in [1], however, we start with simple examples before proceeding to the 
more complex ones. Specifically, in section 2.1, we first consider the simplest 
possible example of f{z) = z and show that it satisfies the generalised root 
identities (1) for arbitrary {zo,n) G C^. In this case, of course, the root side, 
r(sO)M)) is trivial to compute, but the derivative side for arbitrary /j, is already 
interesting, requiring the usual distributional calculations when /x e Z<o and 
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the use of generalised Fresnel integrals when /i G M \ Z. 

We then use this result to extend to the case of arbitrary polynomials in 
section 2.2, and to the case of r(2; + 1) in section 3. 

Finally we consider the case of ( in section 4. By considering the case of 
arbitrary /x and Sq we shall be able to verify numerically that ( satisfies the 
generalised root identities (1) when fj, and sq are in half-planar regions of C, 
and also discuss how these numerical computations may be both extended to 
all of (soiA*) G and used to cast light on certain other theoretical gaps 
identified in [1], in particular the legitimacy of the claim there that we may 
take Clim In (i) = on the root side for ^ for /x = 0, /z = — 1 and fi = —2. 

These results in turn, we believe, give strong numerical evidence in support of 
the idea that ( does satisfy the generalised root identities at arbitrary (and 
in particular at /x = 0, —1, and —2) and hence that the claims in [1] are valid. 

2 The Root Identities for f{z) = z and for Poly- 
nomials 

2.1 The Case of f{z) = z 

In this case / has a single root of multiplicity 1 at 2; = and so the root-side of 
the generalised root identities is given by 

r{zo, IJ,) = -jr (2) 
^0 

Wc now consider the derivative side in several steps. 

(i) The case of /j. G Z>o: When fi = I the derivative side is 

Id 1 

d{zo,l) = WJ\^(^^^)\z=zo = — =r{zo,l) (3) 
1(1) az Zo 

and by direct differentiation it is then trivial to see that we likewise have 

d{zo, /x) = —rr = r{zo, /x) for all /x e Z>o (4) 

so that certainly /(z) = z satisfies the generalised root identities for n G Z^q.^ 

(ii) The case of /x e Z<o: Recalling the definition from [1] of the derivative 

side, that 

1 1 r°° 

and that J^[lnz](^) = — 27r^^^^^ (where HoiO is the odd Hcavisidc function) we 
can calculate the derivative side via the sort of distributional calculations used 
in [1]. For /x = we have 

^Note in passing that these results can be derived via integration from the Fourier definition 
of the derivative side in [1] using Cesaro methods, but there is obviously no need when n 6 Z>o 
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d{zo,0) = A_ 1^ ^ e-o«d^ = 5o(0[e^^°«] = 1 = ^^0,0) (6) 
while for = —1 we have 



d{zo, -1) = ^ 1^ ^ e-««de = z • <5;,(0[e^^°«] = -zo = r{z„ -1) (7) 

and similarly for /i = —2, —3, ... so that /(a;) = z also satisfies the generalised 
root identities for ji S Z<o. 

(iii) The case of fj, G (0, 1): Assume initially that zq is also real and positive. 
Then the derivative side, d{zo,fi) is given by 



2r(M) 



^oCOe^-'CcosC^oO + isin(^oO)d^ 

o 

(1 _ e^'rCM-i)) J^°° ^p-i cos(zoOdC 
+i • (1 + e^'^(''-i)) /~ ^"-1 sin(0oOd^ 



(9) 



on recalling the definition of Ho{^). Changing variables in these integrals by 
letting V = (zo^)^ and noting that dv = iJ,ZQ^^~^d^, it follows that 



e-2f 



But recall the generalised Fresnel integral identity 

'sm{x'')dx=^^^^^^^ (11) 



/ 

Jo 



a 
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defined initially for a > 1 and then extended by analytic continuation.^ 
From this it readily follows also that 

r™(x-)d. = I<i>^5i<M (12) 

Jo 

for a > 1 and extended by analytic continuation. 

Combining (11) and (12) in (10) it follows that we have 



■^When a = 2 this reduces to the classical Presnel integral from optics sin(a;^)da: = \/| 
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dizo,n) = ^J|^{(l + e"-)cos(^) + z.(l-e^-^)sin(^)}.z-'' 

= ^ {(1 + e''^'')(e*5'' + e-'S") + (1 - e''^'')(e'5'' - e-'^'*)} • z^'' 
= —i:r=r{zo,fJ,) (13) 

^0 

Thus we see that f{z) = z also satisfies the generalised root identities for 
arbitrary ^ € (0, 1) and zq € M>o- 

(iv) Extension to arbitrary /i, zq € C: Keeping zq € M>o initially, we can 
then extend from fj, £ (0, 1) to all of /U G M \ Z by a combination of differentiation 
under the integral w.r.t zq and homogeneity arguments (which guarantee there 
are no constants of integration). 

We thus have now verified that f{z) = z satisfies the generalised root iden- 
tities for arbitrary /x e M and Zq G IR>o- By unique analytic continuation they 
therefore continue to hold first for arbitrary n G C and zq € M>o, and finally 
also for arbitrary /x e C and zq e C, except for a branch cut from zq = — oo to 
Zq = when n ^ Z. 

2.2 The Case of Polynomials 

Since f{z) = z satisfies the generalised root identities for arbitrary /i, zq G C it 
follows readily from (5) and elementary properties of the Fourier transform (in 
particular that J^[/(a;-|-a)](^) = ^[/](C)e'"^) that the generalised root identities 
arc also satisfied for arbitrary zo,/j, G C by /(z) = z — a. 

But then it also follows immediately from the fact that \TL((^z — r\){z — ri)) = 
\a.{z — ri) + ln(2; — T'i) and the obvious additivity of the root side that the 
generalised root identities must also be satisfied for arbitrary zq , /i G C by any 
polynomial j?(z). For any given p(z), the structure of the branch cuts in Zq 
when jjb ^ 7j (corresponding to the branch cut on zq G (— oo,0) for f{z) = z 
replicated for each factor) is slightly more complicated, but there are no branch 
cut complications for generic and in particular for all zq in a half-plane to 
the right of the root with the most positive real part, which is all we require in 
certain subsequent arguments. 

3 The Root Identities for T{z + 1) 

In this case r(z + 1) has poles of multiplicity = —1 at —1, —2, —3, . . . and 
so the root side of the generalised root identities is 

.K.„,.) = -e* 2:^^ (14) 
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which is classically convergent for Re{p) > 1 and Cesaro convergent for general 
II (renormalisation being required as per [1] when /x = 1). For example, using 
the Euler-McLaurin sum formula and the same sort of reasoning as used in [1], 
when < /X < 1 we have rr{zo, fi) given by 

rr{zo, n) = -e-'^ lim \ V , ^ - + I (15) 

while for — 1 < < we have 



l-- (^+^)"^ -l(.o + fe)-^;> (16) 



and so on. For ^ Z there is a single branch cut on zq G (— cxd, — 1) (with 
branch points at —1, —2, —3, . . .). 

On the derivative side, using the working from [1], section 3.4, example 2, 
we have 



dr{zo,li) = 



1 



le's 



r(M) 



^o(0 {iO 



H-l 



ei« - 1 
1 - 



2! 3! 



(M) _L Ml! 

2! 3! 



ei(^o+i)€d^ 



(17) 



In [1] we showed already that (ir(^oi m) — ^r(^Oj m) when /i G Z (after renor- 
malisation when /X = 1) and thus that F satisfies the generalised root identities 
when fi is an integer. 

Now, taking the case when S M \ Z, the working in section 2.1, in particular 
(2) and (8) has shown that 



r(p) J- 



(18) 



for arbitrary p € M \ Z and zq ^ (— oo, 0). 

It follows in (17) that, for arbitrary € M \ Z and Zo ^ (— oo, — 1), we have, 
after simplification. 
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1 _ l^iTTfl 1 

' (zo + l)("-i' (^0 + 1) 



_|_iLp*'r(Ai+l) 1 



720 (zo + l)<^' + ^' 



1 _i 1 At-l 
^ + 2 (zo + 1) 



" 2!(zo + l)2 



(M-1)m(m+1)(m+2)B4 
4!(2o + l)^ 



(19) 



>(20) 



as an Euler-McLaurin-style asymptotic expansion for the derivative side. 

Using (14)-(20) we may then write code to numerically check the generalised 
root identities for r(z + 1) when ^ ^ Z. In particular, for any given ^ Z we 
can evaluate drizo^ij) to any desired degree of accuracy for 3f?(zo) sufficiently 
large using sufficient terms in (19); and verifying drizo,fi) = rr{zo,fi) on such 
a half-plane in zq (sufficiently far to the right for the given /z) suffices to verify 
their agreement for all zq by analytic continuation for the given fi. 

Numerical investigations implementing (14)-(20) have been implemented in 
R-code by Dana Pascovici and it has been found that r(2:-|-l) does indeed appear 
to satisfy the generalised root identities for arbitrary zq, /i S C (with branch cut 
on Zq E (— oo, —1) when /i <^ Z). For example, consider the following: 



Re Root and Derivative side Gamma 




muVec 
2000000 roots 
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Im Root and Derivative side Gamma 
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muVec 
2000000 roots 

These pictures show agreement between real and imaginary parts of dr{zQ, fj.) 
and rr{zo,ii) for a range of /i-values from —1 to 4 for zq = 10.381 using 
2,000,000 roots to approximate rr(2o,/^) and truncating the expression for 
dr{zo,^) at ^^^^qr^yx^+sy ■ "^^^ following picture also shows how, for /i small 
(/i — —0.01), the function rr{zQ, —0.01) closely approaches the limiting function 
rr{zQ, 0) which we know from section 3.4 in [1] is given by rr{zQ, 0) = zq + ^. 



Re Root and Derivative sides Gamma as functions of zD 



Re Derlv side 
Re Root side 



~1 1 1 1 1 1 

-0.5 0.0 0.5 1.0 1.5 2.0 2.5 



zOVec 

mu=-0. 01, 2000000 roots 
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The R-script for Gamma is made available with this paper and further nu- 
merical tests may be performed as desired. 



4 The Root Identities for ({s) 

In [1], section 4.1, we derived an expression, using the Euler product formula, 
for the derivative side of the generalised root identities for (, which applies for 
arbitrary so,/i € C, namely 



dd^o,f^)^-^, E (l^P)M \ (21) 



r(M) 

\' / p prime 



For 5f(so) > 1 this is convergent for arbitrary n G C; and for any given jj, can 
then be uniquely analytically continued to arbitrary .so E C from this right half- 
plane in sq. In [1] we confirmed that ^ satisfies the generalised root identities 
when /It e Z>2 (for which the root side is convergent) and that it satisfies them 
for /i = 1 also after removal of an obstruction, with this result for /i = 1 
corresponding to the Hadamard product formula for the related function ^. 
We also noted that (for 5R(so) > 1 initially and then all sq € C by analytic 
continuation) equation (21) clearly implies that 

dciso, fJ-) = as a function of sq whenever /i e Z<o (22) 

In [1], section 4.2, we then performed explicit Cesaro calculations of the root 
side, r{so, /i), for ^ in the three cases fi = 0,-1 and —2. In the first two cases 
we confirmed that 

rc(so,0)=0 and rc(so,-l)=0 (23) 

also, at least modulo an estimate on S{T) in the latter case which holds con- 
ditional on RH; so that C continues to satisfy the generalised root identities in 
these two cases. When = — 2, however, we found that 

r^{so,-2) = ~^^0 (24) 

based on the assumption of the Riemann hypothesis; from which we concluded 
that the Riemann hypothesis appears to be false. At no stage in [1], however, 
did we consider the root side, rQ{so,fj,), for C for arbitrary /j, ^ Z, which we 
noted represented a gap in rigour in our arguments, since care miist b(^ taken in 
ascribing Cesaro limits at isolated values of a parameter (fi) in which analytic 
continuation is being undertaken. 

We thus now consider the possibility of arbitrary /i ^ Z in the root identities 
for In (21) this extension requires no further work for the derivative side, 
which is convergent for arbitrary fj, when K(so) > 1. 
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On the root side, we write (sq, m) ^-^ the sum of three separate contributions 
from T (trivial roots), P (simple pole) and NT (non-trivial roots) as in [1]: 



rdso,IJ')= tt (so ,IJ)+rp (so , /") + tnt (sq , m) (25) 

where 



rp{so,ii) = -e-"^^-^^, and (27) 
rMso,,^) = e^- E TT;^ (28) 

{so-NT} ^ " '^'^ 

When 5ft(/i) > 1, this is also classically convergent for generic sq G C and so 
may be evaluated to any required degree of accuracy by including sufficiently 
many roots, both trivial and non-trivial. Here, by generic sq we mean sq not 
on either (— oo, 1) nor on any of the rays running from a non-trivial root pj = 
A + ill horizontally to the left towards — oo; these being branch cuts of terms 
in the expressions above. Since it is all that is required for our arguments, and 
simultaneously avoids all of these branch-cut issues, from now on we shall always 
restrict consideration to Sq in the half-plane 3?(so) > 1. 



4.1 The Case of 3f?(//) > 1 

For such So we may then readily test numerically whether ( does satisfy the 

generalised root identities also for arbitrary 5ft(/i) > 1 by simply evaluating 
the expressions (21) and (25) to the required degree of accuracy and checking 
whether we do indeed have 



dc{sQ,fJ-) = rt;{so,^J.) (29) 

For the derivative side this requires including sufficiently many primes p (the 
larger sq the fewer primes required to reach a given accuracy); for the root 
side it requires including sufficiently many trivial roots for rx and sufficiently 
many non-trivial roots for tnt (the more positive /U > 1 is, the fewer such roots 
required for a given level of accuracy). 

In order to implement such numerical tests a list of primes and a list of non- 
trivial roots must be used. Such lists can readily be obtained and downloaded - 
e.g. a list of the first 2, 001, 052 NT zeros can be downloaded as a text file from 
[2], and in fact [3] contains files listing up to the first 1, 000, 000, 000 NT zeros. 

Using these we may then readily code up checks for random values of sq > 1 
and /i > 1 and see that ( does indeed seem to satisfy the generalised root 
identities (29) for arbitrary ,so, /i in this region (starting with real sq, H and then 
extending by analytical continuation). 
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Indeed Dana Pascovici's R-script for available with this paper, contains 
such code (currently using the first 10, 000 primes (for d^) and the first 2, 000, 000 
trivial and NT zeros (for r^) but readily adaptable to use more of either if de- 
sired). This thereby allows rapid, systematic checking and on the basis of this 
we can confirm that C does indeed appear to satisfy the generalised root iden- 
tities for arbitrary ^(sq) > 1 when 9fi(/i) > 1, albeit that we cannot extend the 
testing down to /z-values too close to 1 without radically increasing the number 
of roots used in our numerical approximations for r^(so, /x), for obvious conver- 
gence reasons. The following pictures, for example, illustrate the agreement of 
both real and imaginary parts of (i^(so,/^) and r^(so,/i) for sq = 5.1238 for a 
range of /i- values greater than 1.5. 



Re Root and Derivative side Zeta ([nu>1) 



Re Deriv side 
Re Root side 



T 1 1 1 1 r 

1.5 2J} 2.5 3.0 3.5 4.0 



muVec 

2000000 triviai roots, 2000000 NT roots 
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Im Root and Derivative sideZeta (mu>1) 



Im Deriv side 
Im Root side 
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muVec 

2000000 triviai roots. 2000000 NT roots 



4.2 The Case of 3f?(/i) < 1 

When < 1 the evaluation of both tt and tjvt on the root side requires 

Cesaro methods, since the sums in (26) and (28) are no longer classically con- 
vergent . 

4.2.1 Formulae for numerical evaluation of r(^{so, /i) when /i £ (—1, 1)\ 

{0} 

We take the components of the root side in turn: 

Part (a) - rT{so,fi): For rT{so,fi) we use the Euler-McLaurin sum formula 
for the partial sums st- Since the shifted roots occur at sq + 2j, j G Z>i. we 
let z = Sq + 2k + a with < a < 2 and write 



rriso, /i) — e"^ ■ Clim st(soi ^) 



(30) 



where 



(31) 



Now, by Euler-McLaurin, 
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E 



= < 



MM+l)(M+2) (g^^2fc)-^-=^-... 



(32) 



and so we can obtain the Cesaro limit in (30) by the usual Cesaro-style calcula- 
tions on terms of the form (sq + 2k + a)~^~^ (see e.g. [1] or the earlier working 

for r(^ + i)). 

For example, considering just the case < K(/u) < 1 initially, here we have 

(so + 2fc + a)i-^ = (so + 2fc)i-'^|l + (l-/z)^^^^ + ...| 

= {so + 2ky-^ + o{l) (33) 
and thus since pure powers {p ^ 0) have generalised Cesaro limit 0, so 



t't(so,m) = e""'' lim 



l{so + 2k) 



i-p. 



k^oo I ^ (so + 2jy 2 1 - 



(0 < K(m) < 1) 
(34) 

where this is now a classical limit. 

Similar but more involved calculations can be performed on the lower order 
terms to extend (34) to a formula applicable for —1 < 3f?(/x) < 0, then — 2 < 
3?(/x) < — 1 and so on. We shall only consider the next case of — 1 < 5i(/ti) < 
here and in this case the formula becomes 



rT(so,M)=e'"'' 



f 1 

lim < > 

fc^-oo I (so + 



Part (b) - rp(soiA*): This is clearly given by 



r-p(so,M) = -e'"" 



(so - 1)'^ 



(35) 



(36) 



for any Sq and /i. 



Part (c) - rjvx(.so. /'): Here wc use the Riemann-von Mangoldt formula for 
the counting function N{T) and mimic the argumentation in [1]. Write NT = 
NT^ U A^T_ where refers to the subset of non-trivial zeros with positive 

imaginary part and A^r_ to the subset with negative imaginary part. We work 



13 



initially on and then simply state the corresponding formulae for NT- 

which are derived in analogous fashion. For a root Pi = Pi + i^i € NT^ we 
have, on writing A = 5 + e,, that 



1 



(so - (A + ii)y 



(so - ^ - Ei - hi) ^ 



. (so - i - Ci) 



7i 



1 • (so — ^— e.) 
1 — ILl— — * 



2! 

,.jn(j:.+l)(/n+2) {sQ-\-eif 



+ ... 



^ (37) 



and therefore, bearing in mind that any roots off the critical line occur in sym- 
metric pairs either side of it and arguing as in [1], we have 



E 



Mi 



(so - PiY 



-*M(so - \) E{.o-iVT+} ^iTi"""^ 



4S0-JVT. 



2!-^ E{.o-iVT+} ^»7i ^ "Ci 



+1 



. ,L»(j:.+l)(/:t+2)(ao-i)'' 
' 3! 

.M(/^+l)(/:i+2)(so-i) 



3 



We next need to express the partial sums for these series in terms of the 
parameter T using the Riemann-von Mangoldt formula that 



where, as in [1], 



7V(T) = N{T) + ^(T) + -^(T) 

TT 



T T T 1 



(39) 



(40) 



and S'(T) is the argument of the zeta function, and 
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<5(r) 



4 



1 

- tan 
4 



11 7 1 
48T ^ 5760 T3 



+ 



31 1 
80640 T5 



1 

2r 



+ 



T 
2 



{\ M) 



1 

2^5 



(41) 
(42) 



Once we have such expressions we will then re-express them in terms of the 

geometric variable z (as in [1]) and try to evaluate their generalised Cesaro 
limits by removing eigenf unctions z'', z^lnz etc and averaging the remainder 
(inT). 

Note, however, that in combining with N{T) only the sums in (38) that have 
no explicit Ci-dependence will be able to be handled. This is because we certainly 
have no understanding of the functional dependence of any a (T) on T. In part 
for this reason, and also for reasons of numerical tractability, for this paper we 
shall only consider the extension of numerical checks to —1 < 3?(/Lt) < 1, /x ^ Z. 

So now suppose —1 < 3?(^) < 1. Then in (38) we have only the sums 
Y.{so-NT+} ^ili^ and Y.{sQ-NT+} being divergent, the former for all 

such /i, the latter only when K(/i) < 0. Take the case of X]{so-AfT+} ^'^ill^ first. 
Considering just iVr_|_ initially, this generates a contribution up to imaginary 
part T of 



/ t-''dN{t) = T-^N{T) + IJ, [ t-^'-'^N{t) dt 
Jo Jo 



(43) 



+T-^'S{T) + ^ /J t-^'-^S{t) dt ) (44) 

^ +lT-^^5{T) + l^,J^t->^-H{t)dt 

Since 5{T) = O(^), the terms involving 5{T) are classically convergent when 
SR(/u) > — 1. And by arguments analogous to those in [1] the terms involving 
S{T) are also convergent under some suitable power of the Cesaro averaging 
operator P (they are classically convergent when 5i(/i) > and require one 
application of P when — 1 < Ji(/Lt) < 0). 

The only terms where there are explicit eigenfunction divergences that re- 
quire removing are therefore those involving N{T). Here, from (40) we have, on 
writing w = ^ , using integration by parts & changing variables to v = ^ , that 



T-^'N{T) +fi [ t-''-'^N{t) dt = (27r)" 

^0 



[u^ ^Inu 

+^^Io 



+h-' 



' 8 

-1 



u 



dv 
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(l-u)2 \2tv) 



"(WO' 
+C + 0{1) 

Similarly, for X){so-JVT+} 7i~''~^ have contribution 

/ t-^'-'^dN{t) = T-^'-'^N{T) + (/x + 1) / t-^'-^N{t) dt 
Jo Jo 

and the only divergent terms come from the N{T) terms, namely 



(45) 



(46) 



T-^-^N{T) + {ij,+ l) t-''-^N{t)dt =(277)-"-! I " 



+C + o(l) (47) 



It follows from combining these calculations in (38) that, when — 1 < 3?(/Lt) < 
1, the divergent part of J2{so-nt+} (so-Uy '^^ ^^''^^ divgtNT+ (so, T) where 



divgtNT+{so,H;T) = e*>(27r)" 



1 ( T 



+ie^f^(27r)-^-i(so - ^ ^^^1 j^l^^^-C^^ 



(48) 



Since ^ Z it is then clear that divgtNT+ (so, ?") can be expressed exactly 
in terms of powers of the geometric variable z = {sq — \) — iT. In fact 



divgtNT+iso,IJ']z) 



> (49) 



and if /U 7^ in (—1, 1) all these terms have generalised Cesaro limit so that 

ClimdivgtNT.{so, 1^', z) = when /z € (— 1, 1) \ {0} (50) 

2— ^OO 

In other words, for e (—1, 1) \ {0} 
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E 



— = Clim I y 



Mi 



{so — Pi)'^ T^-'oo 1 ^ (sq — pi)l^ 



- divgtNT+ {so, A*; > (51) 



where divgtNT+ {sq, /i; T) is as per (48) and this is in fact a classical Hmit if /x > 
and requires only one application of the averaging operator, Pt, if —1 < fj, < 0. 
In similar fashion, for the NT roots below the real axis, NT_ , we have 



E 



{so-iVT_} 



(so-a)^ 



E— /X 
{so-iVT_} 7i 

+HS0 - i) E{.o-iVT_} 7*"'*"^ 

2! 2^ 2^{so-iVT_} 'i 

/!*-2^2 



2! Z^{so--'VT_} li t 



■ /^(m+1)(m+2) , 
' 3! 



(•So - 2fT,{so-NT_}% 



-fi—3 



X \ \ ' — — 3 2 I 

2J Z^{so-JVT_} 7i -I- 



and thus, by identical reasoning, for /i G (—1, 1) \ {0} the divergent piece of the 
partial sums to "height" T for NT- is given by 



divgtNT.{so,lJ.;f) = e-'5''(27r)-'' < 



1 I if In f 
-ie-'5''(27r)-'^-i(so- |) ' ^ ^ ^ 



.1 fry 



(53) 



For /X ^ this can be expressed exactly in terms of the geometric variable 
z = {so — ^) + iT as 



divgtNT.{so,tJ-;z) = < 



(271 



-(27r)-^{^ 



-(27r)-''-ii(ao-i)(|,) 



(54) 
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with generalised Cesaro limit 0, i.e. 



Clim divgtNT- (so, /i; = when e (— 1, 1) \ {0} (55) 
so that, for /x G (-1,1) \{0} 

V - — , = (7Zim < V - — , - divgtNT. (so, /«; T) > 

{.o6^T_} - ^^)'' u,7r>-f - J 

(56) 

where divgtNT-{so, f^'iT) is as given in (53) and again this is a classical limit if 
/X > and requires one application of the averaging operator, Pj., if — 1 < /x < 0. 

If we write r^vi^ (so- /i) and rjvT_(sO)At) for the expressions on the right in 
(51) and (56) we thus have 

r jvT (so , m) = e*'"'" 7 — s„ = rjvT+ (so , m) + f'iVT. (so , m) (57) 
where 



r-jvT+ (so, /") = e*'''*fiVT+ (so, m) and rjvT_ (so, m) = e"'''fjvT_ (so, m) 

(58) 

and the expressions in (51) and (56) then allow numerical computation of 
rNT{sQ,l^) also for /x e (-1, 1) \ {0}. 

4.2.2 Results of numerical tests for ( for /i e (—1, 1) \ {0} 

Based on the results in parts (a)-(c) in section 4.2.1 we have proceeded to 
implement code for numerical evaluation of (sq, m) and thus for testing of the 
generalised root identities for C for /j, € (—1,1) \ {0}. The following general 
observations apply to these efforts: 

(a) With only 2, 000, 000 NT+ roots, convergence is too slow to obtain rea- 
sonable approximate values on the root side when /i approaches from above. 
This is for the same reason that we couldn't obtain reasonable approximations 
even in the convergent region {fi > 1) when /x approached 1 from above. When 
/X is just above an integer (e.g. 0) the residual piece, even after throwing away 
the divergent terms, converges too slowly (like T^^ or T~'^S{T) when n > 
small) to get reasonable convergence with only 2, 000, 000 NT+ roots. 

However, when /x approaches an integer (e.g. 0) from below, one can get 
accurate numerical results, since on throwing away the extra divergent terms 
that appear when /x passes through the integer value, the remaining residual 
term decays much more rapidly (like T~^~^ or T^^^'^S(T) when /x < small). 

Thus for our numerical tests in this section we consider only /x-regions just 
to the left of /X = 1 and 0. This also helps get accurate evaluations of the 
trivial-roots contributions, rT{sQ,fi), with our choice of 2,000,000 trivial roots. 
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(b) As noted, when < /i < 1 we only need to remove divergent pieces and 
then take classical limits (and in doing this we only subtract those terms in the 
expressions (35), (48) and (53) which are themselves divergent when < /i < 1). 

When — 1 < /I < 0, however, for the NT-root contributions we need to re- 
move the divergent pieces in (48) and (53) and then apply the Cesaro averaging 
operator (in T or T) once to the resulting partial-sum function on the critical 
line, whose step-jumps on crossing roots are now bigger than 1 and grow. Al- 
ternatively this can be viewed as being needed in order to handle the oscillatory 
contributions from the 5'(T)-components of these expressions. 

For fj, small negative and approaching 0, however, we see that the step- 
discontinuities approach height 1 and the resulting partial-sum functions become 
very close to piecewise linear functions with negative slope between root-steps. 
This can be seen clearly, for example, in the following picture: 



This shows the partial-sum function for rjvT_|_ (soi m) when sq — 5.1238 and 
fi = -0.1 on the interval 100,000 < T < 100,005. We see that by using 
the trapezoidal rule between each pair of successive A^T-roots we can readily 
implement code to perform the required Cesaro averaging approximately for 
/z < small, and hence numerically evaluate the contributions from rjvT+ [sq, m) 
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and rNT_ {sq, fi). Using these we are then able to verify whether C does indeed 
continue to satisfy the generahsed root identities also when /i < small and 
whether the root-side of these identities, r^(so, n), does also tend to as /i — > 0^, 
like the derivative side d(^{so,n)- 

(c) With the observations from (a) and (b) in mind, in implementing our 
numerical root-side calculations to test the root-identities for ( we have: 

(i) implemented this within the ^ R-script for < /i < 1 and in order to 
generate pictures of the partial-sum function in T when — 1 < /i < (as for 
example in the figure above), but 

(ii) implemented the full calculations involving Cesaro averaging when — 1 < 
/z < in VBA code in XL2007; this can be adapted to the more powerful 
framework of R, but is sufficient for our numerical experiments in this paper 
and has the virtue of enhanced transparency of the code. 

Both the XL file ("RootIdentitiesZeta_mu_-lToO_TestsB.xlsm") and the 
zeta R-script, along with their required source files, are made available with 
this paper to allow a reader readily to examine the code involved, replicate the 
claimed results and conduct further independent tests if desired (see Appendix 
5.1 for further brief instructions on doing this). 



Results: (a) When < ^ < 1 we get close agreement between d^{sQ,^) and 
(so, fJ.) for fj. > 0.75. For fj. < 0.75, as noted, 2, 000, 000 NT roots is insufficient 
to gain a good approximation to the NT-root contributions to 7'^(so,/i), with 
the situation growing worse as /i — > 0+. The following pictures show how the 
convergence on the root side improves appreciably as /i — > 1^ and agreement 
between d(^(sQ,fi) and r^(so,M) thus becomes progressively better as /i — > 1^: 

Re Root and Derivative side Zeta (0<mu<1) 



Re Deriv side 
Re Root side 



— I — 

0.95 



muVec 

2000000 triuiai roots, 2000000 NT roots 
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Im Root and Derivative side Zeta (0<mu<1) 



Im Deriv side 
Im Root side 



0.75 



0.80 



0.85 



0.90 



0.95 



muVec 

2000000 triviai roots. 2000000 NT roots 



We believe that these provide convincing evidence that ( does indeed satisfy 
the generalised root identities also for < /i < 1; this could be further tested, 
over a wider /i-range, by using more than 2, 000, 000 trivial and iVT-roots. 

(b) When — 1 < /i < we likewise get reasonably close agreement between 
dQ{so,ijL) and rij(so,/i) for /i > —0.2 after Cesaro averaging, but 2,000,000 NT- 
roots is insufficient to obtain reasonable root side values when fi < —0.2: 

Re Root and Deriv Sides Zeta (-l<mu<0) 



9^ 



♦ d_ieta_Re 
■ r_zeta_Re 



-0.22 -0.2 -0.18 -0.16 -0.14 -0.12 -0.1 -0.08 -0.06 -0.04 -0.02 

mu 
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Im Root and Deriv Sides Zeta (-l<mu<0) 



-0.22 -0.2 -0.18 -0.16 -0.14 -0.12 -0.1 -0.08 -0.06 -Ma -0.02 



♦ d_zeta_lm 
■ r_zeta_lnn 



These pictures are again for sq = 5.1238. Here, and in all testing described 
for — 1 < /X < 0, the approximations are obtained using 10,000 primes on the 
derivative side with a truncation threshold of 10 (so that d(^{so, /i) values should 
be very accurate), while on the root side we use 2, 000, 000 trivial roots for rx 
and an average over the next 10, 000 roots starting at the 1, 000, 000*'' A^T-root, 
with 5 subintervals in the trapezoidal rule and an offset of 10^^ (to avoid the 
step-jumps at roots) in the calculation of r^x^ and rjvT_ ■ 

The following tables show further detail on the breakdown of the calculations 
for d(^(so,fj,) and r^(so,/i) when /i is very small, namely /i = —0.001, = 
-0.0001 and fi = -0.00001: (insert tables here) 







5R(rc) 


K(rT) 


5R(r-p) 


5R(rjvT+) 


K(rjvT_ ) 


-0.001 


0.00003300 


-0.00009450 


-3.064361 


-1.001413 


4.348981 


-0.283302 


-0.0001 


0.00000330 


-0.00012270 


-3.062147 


-1.000142 


4.343429 


-0.281263 


-0.00001 


0.00000033 


-0.00012547 


-3.061925 


-1.000014 


4.342870 


-0.281056 







3(rc) 


S(rT) 


S(rp) 


3(r-jvT+) 


3('^JVT_) 


-0.001 


-0.0000001037 


0.0000002974 


0.009627 


0.003146 


737.242 


-737.255 


-0.0001 


-0.0000000010 


0.0000000385 


0.000962 


0.000314 


7,360.36 


-7,360.36 


-0.00001 


-0.0000000000 


0.0000000039 


0.000096 


0.000031 


73,591.4 


-73,591.4 



As /i — > we see that: 

(i) Agreement between d(^(so,fi) and r,^{so,fj,) becomes very close and both 
quantities approach O.'^This provides strong evidence, first, that ( does continue 



•^albeit there seems to be a threshold for the root side reflecting the limits of accuracy we 
can obtain using only 2, 000, 000 trivial roots and an average over only 10, 000 roots starting 
at the 1,000,000"* NT-root in our approximate calculations 
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to satisfy the generalised root identities for /i < and, secondly, that the values 
of (ij(so)At) (and hence also r^{so,fj,)) behave continuously across the integer 
point = 0, with these values approaching the limiting value of when jU = 0, 
as calculated in section 4.1 in [1]. 

(ii) Note in regard to (i) that while the root-side does approach as /j. ^ 0~ , 
further confidence in the correctness of the root identities comes from the fact 
that both the real and imaginary parts of r^(so, /i) arc made up from the addition 
of disparate elements {rr, rp, Tnt^ and tnt-) whose contributions (see the 
above tables) do not themselves become small (and in some cases appear to 
diverge as /i 0~). 

(iii) In this regard note further that as ji goes from —0.001 to —0.00001 
the value of rr more and more closely approaches —\so — \ (as per equation 
(49) for /i = in [1]); the value of rp more and more closely approaches — 1 
(which it equals when /x = 0); and the value of tnt = tnt^ + tnt- more and 
more closely approaches |so + | (as per equation (65) for = in [1]). Thus 
the behaviour of each component of the root side as /i — > 0~ does appear to 
continuously approach the behaviour calculated in [1] for these components in 
the limiting case /i = 0. 

4.2.3 A conjecture regarding the appearance of the lii(|) diver- 
gences in the calculation of rf(so,M) when = and numerical 
test of this conjecture 

In [1] the calculations of the TVT-root contributions to the root sides, r^(so,/i), 
of the C, root identities at /x = 0, — 1 and — 2 all relied on claiming that a value 
of should be ascribed to the 2d Cesaro limit 

Clim In f ^) = (59) 

But in [1] this claim is only justified in a highly heuristic fashion.^ By considering 
the case of non-integer as /x approaches 0, however, we may formulate a precise 
conjecture as to how the expression In (|) arises naturally in the limit and why 
it would have to have generalised 2d Cesaro limit in order to ensure continuity 
of the analytic continuation of ^(^(so, /i) in /i across /i = 0. 

To see this, suppose /x 7^ is very small and consider the expressions 
■rjvx+ (.So, /i) and rjvT_ (so, A*) defined by (51) and (56). With sq real initially 
as usual, it is easy to see from the expression for divgt]^p^{so, ij;T) in (48) and 

for divgtNT- {so, /x; T) in (53) that rjvT+ (so, m) and fjvT_ (so, m) must be complex 
conjugates: 

If fNT+{so,lj) = a-\-ib then fNT-{so,lj) = a — ib (60) 

Since the real part of rNTifyQ, f-i) at fi = is ^.so -|- | (as shown in equation (63) 
in [1]) so, for the real parts, continuity of analytic continuation across /x = is 
equivalent to requiring that 

''or claimed as an implication of C satisfying the /U = root identity, which can then at 
least be legitimately used for the /j, = —1 and ix = —2 calculations 
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lim 3?(e*''^(fjvT+(so,A*) +r'jvT_(so,A*))) = ;^so + ^ 

or, in light of (60) and the fact that lim e'^'^^ = 1, that 

1 3 

^ITO K(fArT+ (.So, M)) = 7«0 + 7 (61) 

As for the imaginary part of r^vT, this automatically tends to since by 
(60) the imaginary parts of rNT+ (so, m) ^-^d tnt^ (sq, m) cancel and e"^** — )• 1 as 
/i ^ as noted. Continuity of analytic continuation across /i = thus requires 
that the overall imaginary part of r^vr also be at /U = 0, and since by equation 
(63) in [1] this is precisely given by 

3(rjvT(so,0)) = -^{so - b ■ Glim In (I) (62) 

so this shows why we must have Clim In (|) =0. 

But we can go further and, by considering tnt^ {so^ m) and tmt^ (so, m) sep- 
arately we can conjecture as to how the In (|) term arises naturally in the limit. 
Recall that (51) may be re-expressed as saying that for /it 7^ small 

7 — j^' y = divgtNT+ (so, M; z) + fNT+ {so, m) + i?(so, m; T) (63) 

where divgtNT+{so, IJ-', z) is as given in (49) and where either R = o(l) (if < 
/z < 1) or P[R]{so,fi;T) ^ as T ^ 00 (if -1 < ^ < 0). 

On the other hand, setting u = ^ and combining result (60) from [1] and the 
further observation from the same section of [1] that P[S]{T) — )• as T ^ 00, 
we have the corresponding claim that at = 

EMi . .X « / Ix, ,1 3, „ „N 

(s -n)o = ^u-lnu-{-+i)u- —{so- -)lnu+{-so + -) + R{so,0;T) 

(64) 

where P[i?](so, 0; T) ^ as T ^ 00. 

But now it is easy to see that, as /z — >■ 0, the first and second terms in (49) 
turn into the corresponding terms in (64) while we saw in (61) that the real part 
of rNT^{so, fi) must limit to the term jSq + |- As for the final term in (49), 
using the the Taylor series expansion 

-u-^' = -{l- n\nu + 0(m^)} = --lnu + 0{n) (65) 
MM M 

we see that, as /x 0, it formally turns into precisely the term —^{sq— 5) Ih'" 
in (64) together with an additional pure imaginary constant term ^{so — 5)^ 
which blows up as ^ 0. 
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In order to match (64) in the hmit as /i — > we would therefore have to have 
?^JVT+ (soj A*) containing an exactly offsetting imaginary component, — 2^(so — 
|)^, and so, combined with (61), we conjecture that 

Conjecture: As 

fNT^{so,^i)-{\so + -^)-^{so-l)^ 

meaning more precisely that 

Wnt+ {so, l^)) = \so + l+ 0(1) (66) 



and 



^{fNT+iso, ^^)) = -^iso-b-+C + 0(1) (67) 

for some constant C. 

In this perspective then, the imaginary parts of r]yT+ (so, m) ^^id r^T^ (so, m) 
are both diverging as fi but in exactly offsetting ways, and the residual 
In (l) divergence which appears in handling the root side of the /x = root 
identity is arising naturally from applying the Taylor series expansion (65) to 
each of the final and terms in the divergent pieces (49) and (54) for 
^jVT+(sOiA*) and rjvT_ (so5M) and then taking the limit as /z — >■ 0. Then In (|) 
naturally has to have 2d Cesaro limit since it arises from a limit in of the 
expression — which clearly has Cesaro limit at all /x.^ 



Results: We tested this conjecture by calculating the values of 3t{rNTjf. {so, At)) 
and 5(fArT+(so,M)) for ^l = -0.001, /z = -0.0001 and ^I = -0.00001 as shown 
in the following tables (for sq = 5.1238 again): (insert tables here) 





K(fjVT+) 


Pred. Leading Order 


-0.001 


2.0328495452 


2.03095 


-0.0001 


2.0310831931 


2.03095 


-0.00001 


2.0309067012 


2.03095 






9(rNT_,.) 


Pred. Leading Order 


-0.001 


745.20948403 


735.90063 


-0.0001 


7,361.1449162 


7,359.0063 


-0.00001 


73,591.493604 


73,590.063 



We clearly see that both parts of the above conjecture appear to be confirmed 
numerically with ^(rNT+{sQ, n)) becoming progressively closer and closer to 

formal L'hopital's calculation, as per section 4.3 of [4], also then recovers jSq + | in 
the limit as /u under this conjecture 
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i*o + I = 2.03095 and '^{rNT+ {sq, m)) diverging as predicted with leading order 
asymptotics more and more closely resembling —^{sq— ^)^- 

As noted this in turn provides a strong justification for ascribing value to 
the 2d Cesaro limit of In (|) when /i = (as done in [1]) in order to retain 
continuity of analytic continuation across /x = 0; and presumably similarly for 
H = —1 and n = —2. 

5 Appendices 

5.1 Tools for performing numerical calculations used in 
this paper 

The numerical results presented in this paper are generated in R-scripts ("R 

code Gamma.R" and "R code Zeta.R") and an XL2007 spreadsheet with VBA 
code ("RootIdentitiesZeta_mu_-lToO_TestsB.xlsm"). The supporting files for 
the zeta R-script are two .txt files ("FirstlOOOOPrimes.txt containing the first 
10, 000 primes and "zero_aII2m.txt" containing the first 2, 000, 000 NT+ roots 
of C as obtained from the zip file contained at [2]). These are also used to 
populate the required areas of the XL spreadsheet in order for its functions to 
work correctly. The spreadsheet has been saved without this data embedded for 
reasons of space, which is also why the file for the 2, 000, 000 iVT_|_ roots is left 
for the reader to unzip and rename as "zero_all2m.txt". For both the R-scripts 
and the XL spreadsheet the code is public and completely transparent. It is 
clear from this code how to run the applications (e.g how to have the functions 
correctly read the input roots and primes) and for the XL spreadsheet there is 
additional brief guidance given within the spreadsheet itself. 

Both applications (the R-scripts and the XL/ VBA) can also easily be adapted 
in order to conduct further testing as desired. 
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